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A search for a manifestly exotic S = —2 baryon state decaying to E~n~ , and its neutral partner 
decaying to S~tt + , has been performed using 220 pb^ 1 of pp collisions at \/s = 1.96 TeV collected by 
the Collider Detector at Fermilab. The S~ trajectories were measured in a silicon tracker before their 
decay, resulting in a sample with low background and excellent position resolution. No evidence 
was found for S = —2 pcntaquark candidates in the invariant mass range of 1600-2100 MeV/c 2 . 
Upper limits on the product of pentaquark production cross section times its branching fraction to 
S~-k + '~ , relative to the cross section of the well established S(1530) resonance, are presented for 
neutral and doubly negative candidates with pr > 2 GeV/c and \y\ < 1 as a function of pentaquark 
mass. At 1862MeV/c 2 , these upper limits for neutral and doubly negative final states were found 
to be 3.2% and 1.7% at the 90% confidence level, respectively. 
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from minimal quark-antiquark or three quark configura- 
tions of standard mesons and baryons, have taken place 
since the outset of the quark model. Until recently these 
searches yielded little evidence for such states. Most no- 
tably in the 1960's and 1970's there were searches for 
a Z*-baryon with positive strangeness in bubble cham- 
ber experiments with incident beams of pions and kaons, 
which ended with negative results. These earlier searches 
are summarized in detailed reviews by Hey and Kelly Q 
and more recently by Jennings and Maltman Q. 

During 2003, however, the situation changed dramat- 
ically. A narrow resonance, decaying to nK + , was re- 
ported by the LEPS Collaboration in their study of 
exclusive photon-nucleon reactions [|[. A state decay- 
ing to nK + must have positive strangeness, an exotic 
quantum number for a baryon, and therefore cannot ex- 
ist in a three-quark model. A pentaquark interpreta- 
tion was employed, suggesting its quark contents to be 
uudds. The mass of the resonance was measured to be 
(1540 ± 10)MeV/c 2 , in remarkable agreement with the 
theoretical prediction for the 8 + state made in the frame- 
work of the chiral soliton model by Diakonov, Pctrov, 
and Polyakov Various experiments, using incident 
beams of real and quasi-real photons, kaons, and neu- 
trons, claimed to confirm the initia l rep orts of + decay- 
ing to nK+ and pK° s @, i 0, |, |, Efl dH, GJ, [3 . 

Another manifestly exotic state, a doubly negative 
baryon decaying to E~tt~ [l4l |. and its neutral partner 
decaying to S~tt + , were reported by the NA49 Collabo- 
ration in y/s = 17.2 GeV pp collisions at the CERN SPS 
[l5| . These states have been tentatively labeled $ and 
<f>° respectively, and the mass of the combined $ __ and 
$° signals was measured to be (1862 ± 2) MeV/c 2 . These 
resonances were interpreted as the I3 = —3/2 and the 
I3 = +1/2 partners of an isospin quadruplet of five-quark 
states with quark contents dsdsu and dsusd, respectively. 
No independent confirmation of these states exists. 

All ostensible pentaquark signals are based on rela- 
tively small data samples. Typically 20 - 100 signal 
events are seen with estimated significances ranging from 
three to five standard deviations. With the exception 
of the ZEUS and HERMES Collaborations, which both 
have reported observations of 9 + in the pKg mode, 
many searches for strange pentaquark states produced 
in high-energy, high-statistics, fixed target or collider ex- 
periments have yielded negative results flrl [l7l. llH , IT9L 

0, m m m, m, m, m m, m, m nonHinHi . Mo re 

recently, a high statistics analysis of the exclusive reac- 
tion 7p — > K a K + n performed by CLAS, which was one 
of the experiments initially reporting the observation of 
the 9 + , yielded no evidence for this state [33j . 

We report the results of searches for exotic baryons 
$ 0,_ ~, decaying to 5~tt + ~, in pp collisions at y/s = 
1.96 TcV. The data were recorded with the CDF II detec- 
tor at the Fermilab Tevatron and correspond to 220pb _1 
of integrated luminosity. 

The CDF II detector is described in detail elsewhere 
[34[. The aspects of the detector most relevant to this 



analysis are the tracking and the three-stage triggering 
systems. The tracking system lies within a superconduct- 
ing solenoid, which provides a uniform axial magnetic 
field of 1.4 T. A system of double-sided micro-strip sili- 
con sensors (SVX II) [35[ is arranged in five concentric 
cylindrical shells with radii between 2.5 and 10.6 cm from 
the beam-line [36| , and approximately 90 cm in length. 
Each silicon layer contains axial strips to measure track 
positions in the plane transverse to the beam-line. Stereo 
measurements of ±1.2° and 90° with respect to the beam- 
line are also available. The remainder of the magnetic 
volume is occupied by a drift chamber, known as the 
central outer tracker (COT) [37|. The COT has a radial 
extent from 40 to 137 cm, and extends 3.1m along the 
direction of the beam. The sense wires are arranged in 
8 super layers of 12 wires each, which alternate between 
axial and ±2° stereo measurements to provide tracking 
in three dimensions. 

The CDF II detector docs not have a dedicated trig- 
ger selection for $ pentaquark production. We therefore 
chose data sets obtained by two complementary types of 
triggers. The first dataset consists of the largest number 
of inelastic pp collisions recorded by the experiment, and 
is based on a displaced-track trigger, even though the i> 
decay products usually produce tracks that have cither 
too small or too large a displacement relative to the pri- 
mary vertex to satisfy the trigger requirements. The sec- 
ond sample is obtained by a calorimctcr-based jet trigger. 
This trigger has the merit of being largely independent 
of the pentaquark production and decay properties rele- 
vant to the displaced-track trigger. However, due to its 
low jet-energy threshold, only a relatively small fraction 
of such triggers could be recorded, thereby limiting its 
statistical sensitivity for this search. 

The displaced-track trigger was designed primarily for 
recording hadronic B-decays. This trigger requires the 
presence in the event of a pair of oppositely charged 
tracks identified in the transverse view by a fast hard- 
ware track fitter called the silicon vertex trigger (SVT) 
[38| . Each track is required to have transverse momen- 
tum pt > 2.0 GeV/c, and an impact parameter (dp) which 
falls within the range 120 fim. < |do| < 1 mm [39J . In ad- 
dition, the intersection point of the triggered track pair is 
required to have a 200 /im transverse displacement with 
respect to the beam-line, the scalar sum of the two trans- 
verse momenta must exceed 5.5 GcV/c, and the azimuthal 
angular difference between the tracks is required to be 
within the range 2° < Acjx 90°. 

The calorimeter based trigger requires the presence of 
a jet with transverse energy Et > 20 GeV (the Jet20 
trigger). The CDF II calorimeter system covers a pscu- 
dorapidity region of 1 77 1 < 3.5 I36| and consists of electro- 
magnetic and hadronic sections [401 ] . The physical towers 
are grouped into trigger towers that span 15° in <f> and 
approximately 0.2 in 77 for central towers (|7j| < 1.0). 
The level- 1 trigger requires a single tower with an Et > 
5 GeV. At level 2, a calorimeter cluster is required to 
have an Et above 15 GeV. The level-3 trigger uses a 
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FIG. 1: Invariant mass spectra from the SVT sample for: (a) An^ candidates that have an associated S~ track in the SVX II, 
where the open histogram shows right-sign A%~ combinations, and the filled histogram shows wrong-sign Air + combinations; (b) 
E~ty + combinations, which display a clear 2(1530) peak; and (c) E~ir~ combinations. Arrows mark the mass at 1862 MeV/c 2 , 
where the NA49 Collaboration reported observing the $. The smooth curves represent fits to the spectra and are described in 
the text. 
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FIG. 2: Invariant mass spectra from the Jet20 sample for: (a) A-k^ candidates that have an associated E~ track, where the 
open (filled) histogram is right-sign Ait~ (wrong-sign Air + ) combinations; (b) E~ir + combinations; and (c) E~tv~ combinations. 
Arrows mark 1862 MeV/c 2 , and the smooth curves represent fits described in the text. 



jet cone algorithm to reconstruct jets with a cone size 
R = \J Ari 2 + Acj) 2 < 0.7 around a seed tower [4l|. Due 
to trigger rate limitations, the jet trigger is prescaled 
such that this dataset contains an effective integrated lu- 
minosity of only about 0.52 pb . 

Pentaquark candidates were reconstructed through the 
decay chain: — > S~ir + -~, E~ — > Att~ and 
A — > ptt~ . Our reconstruction and selection strategy 
was guided by using the ^(1530) — > E~n + signal in 
our data as a model for <I> production. The A candi- 
dates were reconstructed from oppositely charged pairs 
of tracks that had more than 20 axial and 16 stereo COT 
hits and an impact parameter greater than 500 fim. The 
track with the highest transverse momentum in the pair 
was assigned the proton mass. This assignment is always 
correct for true A candidates whose decay products have 
sufficient px to be measured by the tracking system. A 
vertex constrained fit of track pairs was performed for 
the A candidates. The invariant mass was required to be 
within ±5 MeV/c 2 of the world average A mass, and the 
vertex displacement in the transverse plane was required 
to be greater than 1 cm from the beam-line. 

The S~ hyperons were reconstructed in the decay 
chain 3~ — » Att~ by combining A candidates with each 
remaining negative track, assumed to be pion, in the 
event. The track parameters of each three-track combi- 



nation were refit with a E~ decay hypothesis, including 
a mass constraint to the A mass (1115.683 MeV/c 2 ) [42j 
on the tracks originating from the A candidate. Candi- 
dates that were consistent with this hypothesis, and had 
a total invariant mass within 60 MeV/c 2 of the E~ mass 
(1321.31 MeV/c 2 ) [HI, were refit with a E~ mass con- 
straint. The S~ decay vertex was required to be closer 
to the primary vertex than the A decay vertex, and we 
required at least 1 cm radial separation between the two 
decay vertices. 

The lifetime of the S~, ct = (4.91 ± 0.04) cm [Hj], 
and the pr requirements resulting from our acceptance, 
produce a S~ decay vertex that lies outside the inner 
layers of the SVX II for a majority of the events. This 
makes it possible to reconstruct the S~ track from the 
hits left in the layers of the SVX II. The momentum and 
decay vertex of the S~ candidate as reconstructed in the 
COT were used to define its track helix. Hits in the SVX 
II detector made by real S~ hyperons are expected to fall 
near the extrapolation of this helix in the silicon detector. 
A dedicated tracking algorithm was used to attach hits in 
the SVX II to the S~ helix, which required a minimum 
of two axial silicon hits. The hyperon tracks found in 
the SVX II were used, with E~ mass assignment, in the 
subsequent analysis in the same way as standard tracks 
produced by stable particles. 
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The Att~ invariant mass spectrum of £~ candidates 
found in the SVT dataset is shown in Fig. QJa). A £~ 
signal containing about 36,000 candidates is evident. The 
invariant mass spectrum of wrong-charge combinations, 
Att + , shown by the filled histogram, provides a good de- 
scription of the right-sign combinatorial background. No 
artificial signal was created by the selection criteria or 
the hyperon tracking procedure. 

To search for new states, each £~ candidate was com- 
bined with all remaining tracks. These tracks were as- 
signed the pion mass and were required to have pr > 
0.4 GeV/c, at least 20 axial and at least 20 stereo COT 
hits, and 3 or more axial SVX II hits. In addition, £~ 
tracks and pion tracks were required to have impact pa- 
rameters that were smaller than 0.015 cm and 0.1 cm re- 
spectively. The remaining £~7r ± pairs were subjected to 
a common vertex fit with the total momentum of the pair 
constrained to point back to the primary vertex in the 
transverse plane. A vertex fit quality requirement was 
imposed, so that only well measured candidates were re- 
tained for further analysis. 

No additional kinematic cuts were applied to £~7r ± 
pairs. As described earlier, such states are unlikely 
to produce track pairs satisfying the SVT trigger cri- 
teria. Therefore, £ _ 7r ± candidates were not required 
to satisfy the SVT trigger. Since the SVT trigger se- 
lects events containing particles produced in the decays 
of long-lived and b hadrons, our source of pentaquark 
candidates would be predominantly the fragmentation of 
the available light quarks in the event, as is the case for 
£(1530) baryons. 

The invariant mass spectra of £~7r + and £~7r~ com- 
binations are shown in Figs. HJb) and (c) respectively. A 
prominent peak corresponding to the decay £"(1530) — > 
£~7r + is clearly visible in the £ _ 7r + mass distribu- 
tion. No other narrow structures can be seen in cither 
distribution. The smooth curves superimposed on the 
£~7r ± spectra show the results of unbinned maximum- 
likelihood fits to the data. The fit function included a 
Breit-Wigner resonance shape convoluted with a Gaus- 
sian resolution function describing the £"(1530), plus a 
third degree polynomial multiplied by an exponential as 
a model for the background. The width of the Gaus- 
sian was allowed to float, and the intrinsic width of the 
£(1530) was Gaussian-constrained in the fit to the world 
average value of (9.1 ± 0.5) MeV/c 2 0. The fit yielded 
1923 ± 80(stat) candidates above the background at a 
mass of (1531.7 ± 0.3(stat)) MeV/c 2 , in excellent agree- 
ment with world average mass of the £(1530). The fitted 
mass resolution of (3.5±0.6(stat)) MeV/c 2 is in agreement 
with the expected mass resolution based on the full de- 
tector simulation. No excess of events appears at a mass 
of 1862 MeV/c 2 , the mass claimed by NA49 for the $ 
states. Since the production mechanism of pentaquark 
states is unknown, it is possible that the SVT trigger 
bias may preferentially reject events with pentaquarks. 
To address this issue, an identical analysis has been per- 
formed using the data obtained with the Jet20 trigger. 



The invariant mass spectra of Air~ and Att + combi- 
nations reconstructed in the Jet20 data and having an 
associated SVX II hyperon track are shown in Fig. [21a). 
The jet-trigger data sample is significantly smaller than 
the data sample obtained with the trigger on displaced 
tracks. However, this sample of 4870 ± 128(stat) £~ 
candidates is still a factor of 2.5 times larger than the 
sample used by the NA49 experiment in their analysis. 
The £~7T + and invariant mass spectra are shown 

in Figs. [5fb) and (c). Again, the baryon £(1530) is 
seen, and no narrow structures are visible at the mass of 
1862 MeV/c 2 . The £ _ 7r + invariant mass spectrum was fit 
in a similar fashion as previously described with the ex- 
ception that the Gaussian resolution of the £(1530) peak 
was constrained to (3.5 ± 0.6) MeV/c 2 . The fit yielded a 
signal of 313 ± 28(stat) £(1530) candidates. 

There is no indication of S = — 2 pentaquark produc- 
tion in these data. We determined an upper limit on the 
4> pentaquark cross section times the branching fraction 
K($ — > £ - 7r ± ) relative to the £(1530) cross section in pp 
collisions. Knowledge of this ratio can be an important 
input for pentaquark production models. Since the two 
states share the same final state, many sources of system- 
atic uncertainty cancel out in this ratio. We measured 
this ratio for the $ mass of (1862 ± 2) MeV/c 2 , the value 
reported by the NA49 Collaboration, and performed a 
scan of this ratio in the mass range [1600-2100] MeV/c 2 . 

Relative acceptances and reconstruction efficiencies for 
£ _ 7T ± final states of different invariant mass have been 
studied with pp — > 66£(1530) X events produced with 
the pythia event generator [43| and subjected to the 
full GEANT simulation of the CDF II detector, trigger, 
and reconstruction. The pentaquark states were repre- 
sented in the PYTHIA decay table as £(1530) with various 
masses. The pentaquarks were not required to originate 
from hadrons containing a & or quark. The simulation of 
the associated production with bb pairs is a purely prac- 
tical approach to satisfy SVT trigger conditions. Four 
samples were generated for masses of 1530, 1696, 1862, 
and 2027MeV/c 2 in order to obtain the mass dependence 
of the detector resolution and detector efficiency for the 
whole mass range of interest. The generated pr spectrum 
of the £(1530) was re-weighted to reproduce the pr dis- 
tribution of the £(1530) measured in the SVT trigger 
data [Hj]. The latter is shown in Fig. 03 The smooth 
curve shows a fit with an exponential convoluted with 
a Gaussian. The same weight function was applied to 
the pt spectra of the generated pentaquarks in all four 
samples. Lacking any established model for pentaquark 
production, we assumed that the px spectrum of the 4> 
pentaquark is the same as that of the £(1530), other than 
the kinematic effects of heavier masses in the pythia 
generation. 

The sensitivity of the experiment to £~7r ± final states 
is limited by our tracking acceptance to particles that 
have transverse momentum greater than 2 GcV/ c and 
\y\ < 1- Compared to the £(1530), a pentaquark with 
a mass of 1862McV/c 2 would have a factor of 1.5 higher 
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SVT trigger 


Jet20 trigger 


N(S-) 

JV (15.30)) 


35722 ± 326 

i r\o o i on 

1923 ± 80 


4870 ± 122 

o i o _i_ no 

313 ± 28 


a(pp-><I. u X).B($ u ^H-^) 

S(1530)X) 1 J 
_ $— X ) • B($— - S-7T-) 

S(1530)Z) 11 


< 3.2 (5.8) 

< 1.7 (3.1) 


< 3.0 (9.2) 

< 3.2 (10.1) 



TABLE I: Event yields of 5~, S(1530), and upper limits on a(pp -> ^° ~~X) ■ 2?($° •"" -> S~n + '-)/a{pp -> S(1530)X) for 
m($) = (1862 ± 2) MeV/c 2 , at 90% confidence levels obtained by assuming zero natural width of the pentaquarks. Numbers in 
parentheses were calculated with the assumption of a 17MeV/c 2 natural width. 
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FIG. 3: Reconstructed transverse momentum spectrum of the 
£"(1530). The fit function is described in the text. 



detection efficiency at pt = 2 GcV/c; this ratio of effi- 
ciencies decreases with increasing pr. We estimate that 
the ratio of detector efficiencies averaged over the range 
of pr within our acceptance is 1.29 ± 0.05(stat), where 
the uncertainty is due to the size of the simulated event 
sample. There is an additional 10% systematic uncer- 
tainty due to the parameterization of the S*(1530) trans- 
verse momentum. This value is obtained by varying the 
parameters of the fit to the observed S"(1530) pr spec- 
trum within statistical errors, and re- weighting the input 
spectrum of the simulation accordingly. 

The limits on a(pp -> $ ' — X) ■ -> 
S~n + '~)/<r(pp — > S(lb3Q)X) were obtained by intro- 
ducing an additional signal term for the hypothetical 
pentaquark signal, represented by a Breit-Wigner res- 
onance shape convoluted with a Gaussian resolution, 
into the previously described fit function. The ex- 
pected detector resolution at ■m(E~n ± ) = 1862MeV/c 2 
is (6.6±0.7(stat)±0.5(syst)) MeV/c 2 , which was obtained 
through detector simulation. The presence of a $ signal 
was evaluated for two hypotheses for its natural width: 
a natural width of zero, and a natural width equal to 
17McV/c 2 , the latter value being the measured full width 
at half maximum of the peak reported by the NA49 Col- 
laboration [ll|. The values of the relative efficiency and 
the mass resolution at m($) = 1862MeV/c 2 were allowed 
to vary by their uncertainties by adding a Gaussian con- 
straint to the likelihood function. The widths of the 
Gaussians were set to the sum in quadrature of the sta- 



tistical and systematic uncertainties. We determined the 
limit on the production of the $ observed by the NA49 
Collaboration by imposing their measured mass with a 
Gaussian constraint of 2McV/c 2 in the likelihood. The 
limits are established using a Bayesian approach, with a 
uniform prior for the relative production cross section, 
and taking into account S(S"(1530) — > S~n + ) = 0.67 
[42l ]. The limits at the 90% confidence level are summa- 
rized in Table HI 

Available predictions for S = — 2 pentaquark masses 
vary over a range from 1750 to 2027 MeV/c 2 0,111, 
Since we failed to sec evidence for a signal at the mass 
stated by the NA49 Collaboration, we have also derived 
a limit on the relative pentaquark production rate as a 
function of mass in the range [1600-2100] MeV/c 2 in order 
to provide a generic bound on $ pentaquark production 
models. The procedure to determine the upper limit on 
the relative cross section for a pentaquark mass within 
the range [1600-2100] MeV/c 2 was similar to the one used 
for testing the NA49 Collaboration's observation. We 
performed a fit for the relative cross section over this 
mass range in 2MeV/c 2 steps in pentaquark mass. The 
relative detection efficiency and mass resolution entered 
as a function of mass into the fit function. The results 
of the relative cross section scans for neutral and doubly 
negative pentaquarks are shown in Figs. 0] and [SI Neither 
the SVT trigger data nor the Jet20 trigger data show 
evidence for pentaquarks over this broad mass range. 

In conclusion, we have performed a search for exotic 
S = —2 baryon states decaying to 5~tt~ and S~tt + over 
the mass range of [1600-2100] MeV/c 2 . No evidence for 
a signal has been found, and we have set upper limits 
on the relative cross section times branching fraction of 
exotic cascade baryons with respect to the well estab- 
lished resonance £"(1530) for px > 2GcV/c and \y\ < 1. 
We summarize our results by quoting the higher statis- 
tics analysis with the displaced-track trigger, in which 



a(pp -► $X) • 6($ u 



■)/<r(pp->E(1530)X) 



for neutral and doubly negative states were respectively 
found to be less than 3.2% and 1.7% at the 90% confi- 
dence level, for the mass reported by the NA49 Collab- 
oration [l5j]. Consequently, we find that these data are 
unable to confirm the existence of exotic S = —2 baryons. 
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FIG. 4: Upper limits from SVT data (90% confidence level) for the production cross section of $ pentaquarks times the 
<E> — > Sir branching ratios relative to the S(1530) cross section as a function of: (a) <3>° mass, and (b) mass. The solid 

line shows the results obtained assuming zero natural width of the pentaquark, and the dashed line shows the results obtained 
assuming T = 17 MeV/c 2 . The arrows mark the position of the NA49 peak. 




FIG. 5: Upper limits from Jet20 data (90% confidence level) for the production cross section of 4> pentaquarks times the 
$ — > Sir branching ratios relative to to the H(1530) crossection as a function of: (a) $° mass, and (b) $ mass. The solid 
(dashed) line shows the results obtained assuming r = (17) MeV/c 2 . The arrows mark the position of the NA49 peak. 
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